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Abstract 


Stanene is a quantum spin hall insulator and a favourable material for electronic and 
optoelectronic devices. Density functional theory (DFT) calculations are performed to study 
the band gap opening in stanene by investigating the effect of beryllium and magnesium doped 
stanene single layer to study the electronic and structural properties in stanene. The electronic 
band energy of pure stanene without spin orbit coupling (SOC) appear to show no energy gap 
at the Fermi level showing that stanene is a gapless material with Dirac cones at the K point 
and the band gap opens by a gap of 0.08 eV is opened at the K point. The electronic structure 
of Be and Mg doped stanene shows that the Fermi level is shifted towards the valance band 
edge when compared to pure stanene. We have considered 6.25, 12.5, 18.75 and 25% of both 
Be and mg doping. The electronic structure of Be doped stanene show that the Fermi level is 
shifted towards the valance band edge when compared to pure stanene. The Dirac point of 
stanene locates at I’ shifted by 0.38 and 0.5leV for 6.25 and 12.5 %, an energy band gap of 
0.27 and 0.50 eV were obtained above the Fermi level for 6.25 and 12.5% respectively. In the 
case of Mg, the bandgap remains slightly above the Fermi-level and amounts to 0.34 eV for 
6.25 % and 0.43eV for 12.5 %, in the case of 18.75 and 25 % the Dirac cone disappear 
completely, an energy gap of 0.28 eV and 0.60 eV were obtained above the Fermi level for 
6.25 and 12.5% respectively, our findings show that the band gap of stanene open at 12.5% 
doping concentration of both Be and Mg impurities. These obtained band-gap value seem to 
be sufficient for use of alkaline earth metal doped stanene in optoelectronic and such 
applications where stanene is incapacitate for its use to switch on/off devices. 


Keywords: Stanene; doped; Density functional theory; generalized gradient approximation; Quantum ESPRESSO 


VOLUME 01, ISSUE 01, 2021 1 ©DOP_KASU Publishing 


PHYSICSAccess 
ISSN Online: 2756-3898, ISSN Print: 2714-500X 


I. INTRODUCTION 


H ollowing the recent technological advancements 
produced by graphene because of its exceptional 
electronic and optical properties, there has been ongoing 
increasing interest in the research of two dimensional (2D) 
nanomaterials [1, 2]. After, the synthesis of graphene, the 
heavy group-IV elements (Si, Ge, and Sn) are receiving 
considerable attention due to their similar electronic properties 
to that of graphene and broad application prospects that range 
from next generation energy efficient Nano-electronic devices 
[3 4 and 5]. For example, linear band dispersion at the Fermi 
energy, which leads to the charge carriers, behaves like 
massless Dirac Fermions [5,16 and 29]. These massless Dirac 
Fermions propagate with a very good Fermi velocity and show 
high carrier Mobility’s at room temperature. 

Nowadays, a new 2D group IV material, stanene, which is 
a buckled honeycomb structure of a 2D hexagonal tin film, has 
incited noteworthy interest due to its intriguing prospect as a 
quantum hall insulator, topological insulator and topological 
superconductor [6-8]. Stanene a bi-atomic layer of a-Sn (111), 
which has a graphene-like honeycomb structure [9-10]. 
However, the structure is most stable with a_ buckled 
configuration due to weak a-2 bonding between Sn atoms 
[11]. Theoretical studies demonstrated that stanene has a zero- 
band gap without inclusion of spin orbit coupling (SOC), 
whereas it shows a band gap of 0.1 eV after inclusion of SOC 
[12-13]. Stanene behaves like a quantum spin hall insulator 
(QSHI) and its derivatives show properties of large-gap 
quantum spin Hall insulators (QSHIs) [12]. Thus, such QSHIs 
show dissipation less conductor at room _ temperature. 
Furthermore, such QSHIs show excellent thermoelectric and 
quantum anomalous Hall (QAH) effect near room temperature 
[14]. On the other hand, doped QSHI materials show time- 
reversal-invariant topological superconductivity. Thus, doped 
stanene could be a promising material for various applications 
[11]. On the other hand, the zero-band gap in stanene also 
limits its applications towards semiconductor based 
microelectronic devices [15-16]. Therefore, the expectation of 
a band gap opening in stanene has spurred intense scientific 
interests. Various approaches have been taken to modulate the 
electronic properties of stanene. Furthermore, the substrate 
plays an important role to tune the electronic properties of 
stanene due to the strong interaction [17 and 18]. Electronic 
properties of stanene are affected by the substrate because it 
induces band inversion in the bonding and antibonding states 
of stanene, which opens up the band gap in stanine [19]. 
Furthermore, van der Waals heterostructure of stanene 
modulates its electronic properties [20]. Recently, Khan et al. 
[21], demonstrated a heterostructure of stanene with h-BN 
sheet induces a tunable gap in stanine [21]. On the other hand, 
doping is also an efficient way to tune the electronic properties 
of 2D nanosheets by breaking the symmetry of the sublattice 
[11-12, 15]. In this work, we have performed a First-Principles 
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calculations of stanine doped group IV elements (Be and Mg) 
to tune the electronic properties of stanene manuscript. 


Il. COMPUTATIONAL DETAILS 


Our calculations were carried out within the density 
functional theory formalism (DFT) [22-23], as implemented 
in the Quantum ESPRESSO code [24]. We used the 
generalized gradient approximation (GGA) with the Perdew- 
Burke-Ernzerhof (PBE) form to treat the exchange-correlation 
functional [25]. The norm-conserving Trouiller-Martins 
pseudopotential. In order to obtain the optimized 
configuration, the unit cells of considered Single layer were 
relaxed in axial direction until the stress components were less 
than | GPa. The conjugated-gradient optimization of atomic 
positions was allowed to proceed without any symmetry 
constraints until the forces acting on each atom were less than 
0.02 eV/A. Also, the convergence criteria for total energy 
calculations was set at 10-4 eV. We used 50 k-points 
sampling according to Monkhorst-Pack approach for 
calculation of the total energy and band structure [26]. Also, 
XCrySDen (visualization for electronic and _ structural 
analysis) program was employed for visualization of the 
volumetric data such as electron/nuclear densities [27]. The 
Brillouin zone in our calculations is integrated with a 5 x 5 x 
1 and 15 x 15 x 1 Gamma-pack k-point grids for geometry 
optimization and density of state calculations, respectively. 
The vacuum space of 20 A is used in the z-direction to avoid 
any interaction between the adjacent repeating layers. Stanene 
possesses a low-buckled honeycomb structure, and this non- 
planar geometry 1s the most stable structure. Stanene consists 
of tin atoms arranged in a single, hexagonal layer. The buckled 
geometry of stanene leads to the weakening of the m — a 
bonding between tin atoms. This buckling and consequently 
bond weakening increase the overlap between z and o orbitals, 
leading to the increase in the stability of the structure. 4 x 4 x 
1 supercell of stanene was used for the electronic structure 
calculations. Fig. 1 displays the optimized structure of the 
pristine stanene supercell in different views. Beryllium (Be) 
and Magnesium (Mg) single atoms were substituted into the 
stanene single layer to model the doped stanene structures, and 
the effects of theses dopants on the structural and electronic 
properties were examined in detail. During the optimization of 
the structure, all atoms of the system were allowed to relax. 


Il. RESULT AND DISCUSSION 
A. Structural Stability 


Number The stability of stanene single layer usually 
achieved from its stable buckled geometry when compared 
with that of grapheme and germanene, this is because Stanene 
normally has a low-buckled non-planar honeycomb structure, 
and this non-planar geometry is the most stable structure. 
Stanene consists of tin atoms arranged in a single, hexagonal 
layer. The buckled geometry of stanene leads to the weakening 
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of the a — a bonding between tin atoms. This buckling and 
consequently bond weakening increases the overlap between 
m and o orbitals, leading to the increase in the stability of the 
structure. In this work an initial four atoms hexagonal unit cell 
of statnene with space group P63 /mmc (194) and crystal 
parameters of a = b = 462 A, c = 6.68 A has been increases 
and a 4x4x1 supercell of 64 atoms was created by means of 
XCrysden software as shown in Fig. 1. 


Buckling size 


°%o b¢ 04 bel 


Fig. 1(a) The optimized 4x41 supercell of stanene layer and (b) the 
buckling size. 





After optimization by means of relaxation calculation 
within BFGS as implemented in Quantum ESPRESSO, we 
found that the optimized unit cell parameters and Sn-Sn bond 
length are a = b = 4.67 A and 2.63 A, respectively, this 
increase of 0.05 in the parameter may attribute due to the 
effect of spin orbit coupling (SOC) in the DFT pseudopontials. 
This result is in agreements with the number of DFT report 
[11-13]. 

In this work we adopt the model of substitutional doping 
along the reference point (H, M and K) shown in Fig. 1 (a), 
our reference point was determined by chosen the one final 
energy as the most stable position. Table I and Fig. 2 shows 
the calculated total energy of the chosen points and optimized 
doped structures. 


TABLE I CHOSEN DOPING REFERENCE POINT AND CALCULATED 
TOTAL ENERGY: 

4x4x1 supercell Fermi energy (eV) 

of stanene 


Total Energy (Ry) 


Along H-site -2655.35966414 4.2178 
Along-M site -2948.12311484 4.3210 
Along-K site -2811.25381884 4.3379 





Fig. 2 (a) optimized Be-doped stanene (b) optimized Mg-doped 
stanene, all at reference point. 
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From Table I, H-site has been the chosen reference point for 
both the dopants (Be and Mg). First of all, the formations 
energies with the two dopants (Be and Mg) has been 
calculated in order to investigate which among them will be 
more stable using (1). 
Ep = Epure — (NsnUsn + Udopant) (1) 

where Ep is the formation energy, Epure iS the total 
calculated energy of the pure stanene system ns, the number 
of Sn atoms in the unit cell, Ws, 1s the chemical potential of 
Sn atom which is defined as energy per Sn atom in a perfect 
stanene and Ugopant 18 the chemical potential of the dopants 
(Be and Mg) according to the energy per Be and Mg atom in 
Hexagonal Be and Mg. 


TABLE II CALCULATED FORMATION ENERGIES OF THE DOPED 


SYSTEMS: 
Doped systems Formation Energy (eV) 
Be-Doped Stanene 1.318 
Mg-Doped stanene 1.425 


From Table II, the Mg-doped system has higher value than 
that of the Be-doped system, indicating that Be-doped system 
will be more stable than that of the Mg-doped, and this 
increase of positive formation energy may attribute due to the 
different in atomic radius of the dopants. 

The doping concentrations has been investigated by varying 
the dopant as substituting one atom of Sn by one atom of the 
dopants (Be and Mg) as 6.25%, two atoms of Sn by two atoms 
of the dopants (Be and Mg) as 12.5%, three Sn atoms by three 
atoms of the dopants (Be and Mg) as 18.75%, and four Sn 
atoms by four atoms of the dopants (Be and Mg) as 25%, To 
investigate the structural distortion in the differently doped 
systems, the average bond lengths and total energy of the pure 
stanene and that of the differently doped systems were 
calculated and summarized in Table III. 


TABLE III CALCULATED BOND LENGTHS OF THE DOPED SYSTEMS: 


Percentage Bond length Sn-Sn(A) — Bond Bond 

of Doped length Sn- length Sn- 

systems Be-Doped Mg-doped Be (A) Mg (A) 
6.25% pF lls | 2.5382 1.9500 2.2300 
12.5% 2.5021 2.5201 1.9959 2.4655 
18.75% 2.5020 2.5200 1.9999 2.4768 
25% 22D 24112 2.1039 2.6331 


From Table III, a Single Be and Mg (6.25%) doping at Sn 
reference site induces decreases in the Sn-Sn bond lengths 
compared with that of the pure stanene single layer in both 
dopants with values 2.5151 A and 2.5382 A, while Sn—Be and 
Sn-Mg bond lengths become 1.95A and 2.23A, respectively. 
This decrease in Sn-Sn bond lengths may be due to the 
difference in ionic radius between Be and Mg, so substitution 
of Be and Mg at Sn site induces relatively small local structure 
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modifications of such 2D materials this results 1s in consistent 
with experimental results [30] 

On increasing the Be and Mg doping concentration to 
12.5%, Table III shows that the Sn-Sn is 2.5021A and 
2.5201A, for the case of Be-doped and Mg-doped respectively 
showing more decreases as doping concentration increases, 
while Sn—Be and Sn-Mg bond length becomes 1.9959 A and 
2.4655 A respectively, being greatly enlarged compared with 
the corresponding bond length of 6.25% doped stanene 
system. having Be and Mg doping concentration increases to 
18.75%, showed nearly the same bond length as that of 12.5%, 
however the Sn—Be and Sn-Mg bond lengths was also slightly 
increases to 1.9999 A® and 2.4768 A°. When four Be and Mg 
atoms were induced at the Sn sites of the stanene supercell, the 
doping concentration becomes 25% and the Sn-Sn bond 
lengths became 2.225 A and 2.4112 A (showing the greatest 
shrinkage among all the doped system) while the Sn—Be, Sn- 
Mg bond lengths became 2.1039 A and 2.6331 A showing a 
huge increase that making the layer of Stanene so distort. 


B. Electronic Properties 


The calculated electronic properties of pure stanene and 
stanene doped Be and Mg are computed within PBE 
approximation based on DFT along the high symmetry point 
of the first Brillouin zone (BZ). The calculated band structure 
of pure stanene are presented in Fig. 3 (a) and (b) show the 
graphs of band structure of stanene without and with spin orbit 
coupling (SOC). The electronic band energy of pure stanene 
without SOC appear to illustrate no energy gap at the Fermi 
level showing that stanene is a gapless material with Dirac 
cones at the K point and these findings are in good agreement 
with previous results [11, 28]. Conversely, the band gap opens 
by including spin orbit coupling (SOC). Once SOC is 
considered, a gap of 0.08 eV is opened at the K point, which 
is in good agreement with previous results [11,13]. 


Energy (eV) 
Energy (eV) 





a M K OS M K r 


Fig. 3 Calculated Bands Energy of pure stanene: (a) without SOC (b) 
with SOC 


The results of total density of states (DOS) and partial density 
of states (PDOS) of pure stanene are presented in Fig. 4 (a) 
and (b). The results of DOS help to further elaborate the nature 
of the band gap while PDOS gives information about the 
origin of bands. The orbital character reveals that the valence 
band maximum (VBM) and conduction band minimum 
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(CBM) are formed by p- and s-orbitals of Sn atoms. 
Furthermore, the s-orbital of Sn atoms contribute the most to 
the conduction bands than the p-orbital. Also, s-orbital of Sn 
atoms are predominantly to the valence bands near the Fermi 
level. Fig. 4 (b) shows that d-orbital of Sn has no contribution 
for both conduction bands and valence bands. Our results are 
in good agreement with previous first principles investigation 
[1, 8 and 11]. 


(a) Stanene without SOC (b) 
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Fig. 4 (a) The calculated density of states of pure stanene (b) The 
calculated partial density of state (PDOS) of pure stanene. 


C. Doping 

One of the effective ways of tuning electronic properties of 
materials is substitutional elemental doping. Being stanene 
has similar electronic properties as graphene; band gap 
opening in stanene could be possible by _ elemental 
substitutional doping. The two atom of primitive unit cell of 
stanene is extended to 4 x 4 x 1 supercell containing 64 atoms. 
Sn is a group IV element, and if we dope a Be and Mg then it 
induces a hole/electron in the system. Thus, it will break the 
symmetry of the sub lattice, which in turn, might change the 
electronic properties of the sheet. 


D. Stanene doped Be 


Beryllium (Be) doped systems are widely investigated and 
found to be promising for electronic and related applications. 
The previous report demonstrated that introducing Be to 
graphene opens up its band gap, indicating that graphene 
doped Be is a semiconductor material with a band gap of 0.298 
eV [22]. Here, we have considered 6.25%, 12.5%, 18.75% and 
25% of Be doping i.e. one Be atom, two Be atoms three Be 
and four Be atoms is doped in the 64 atomic unit cell of 
stanene. Since higher doping concentration of Be modifies the 
structure of graphene [3], higher concentration of Be is also 
considered in order to modify the structure of stanene. The 
electronic band structure of Be doped stanene (see Fig. 5 (a) — 
(d)) show that the Fermi level is shifted towards the valance 
band edge when compared to pure stanene and this confirms 
that the Be doped stanene is p-type of material. The Dirac 
point of stanene locates at I shifted by 0.38 and 0.5leV for 
6.25 and 12.5 % respectively as shown in Fig. 5 (a) and (b). 
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Interestingly, an energy band gaps of 0.27 and 0.50 eV were 
obtained above the Fermi level for 6.25 and 12.5% 
respectively, which indicate that the Be doped stanene could 
be a degenerate semiconductor. Furthermore, in the case of 
18.75 % and 25% the Dirac cone vanished. Therefore, our 
findings show that the best doping concentration of Be is 
12.5%, because large energy gap was obtained at that 
percentage. 


Gt 
NS 


MN \g 





a 
i 


Fig. 5 Calculated Bands Energy of Stanene doped Be: (a) 6.25% (b) 
12.5% (c) 18.75% and (d) 25% 


The DOS (Density of state) and PDOS (partial density of 
states) for the Sn Be atoms of the Stanene doped Be are 
investigated through the GGA as illustrated in Fig. 6 (a) and 
(b). The DOS and PDOS plots in Fig. 5(b) clearly display that 
there is an energy range (below the Fermi, EF as well as above 
it). Furthermore, the DOS and PDOS plots in Fig. 6 reveal two 
main regions one in the VB and the other in the CB. The first 
region below Fermi consists of the contributions by the s- and 
p-orbitals of Sn atoms and s-orbital of Be atoms. The 
conduction band region is mainly contributed by the p-orbitals 
of both Sn and Be atoms. It 1s also observed that the 
contribution of p-orbital of Be is more than that of Sn to the 
conduction band while the s-orbital of Sn atoms are 
responsible for core bands. 
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(a) Stanene doped Be (b) Stanene doped Be 


Total Density of State 
Partial Density of State 





é 4 -2 0 5 , O 
Eneroy (eV) Energy (eV) 


Fig. 6 Calculated DOS (Density of state) and (b) PDOS (partial 
density of states) for the Be atoms doped Stanene 


E. Stanene Doped Mg 


The electronic properties of Stanene doped Mg in the 
hexagonal structure were executed based on PBE 
approximation. The band structure calculations are performed 
along special line of high symmetry points in the irreducible 
Brillouin zone T (0 0 0), M (1/3 0 0) and K (1/2 0 0) k-space 
by setting Fermi level at 0 eV on energy scale as shown in Fig. 
6, where the red dash line represents the Fermi level. The band 
structure of Stanene doped Mg was also computed for 6.25%, 
12.5%, 18.75% and 25% of Mg doping i.e. one Mg atom, two 
Mg atoms, three Mg atoms and four Mg atoms is doped in the 
64 atomic supercell cells of stanene. The bandgap remains 
slightly above the Fermi-level and amounts to 0.34 eV for 6.25 
% and 0.43eV for 12.5 %, in the case of 18.75 and 25.0 % the 
Dirac cone disappear completely. Conversely, an energy gap 
of 0.28 eV and 0.60 eV were obtained above the Fermi level 
for 6.25 and 12.5% respectively, showing that Mg doped 
stanene could be a degenerate semiconductor. Also, 1n the case 
of 18.75 % and 25 %, the Dirac cone vanished. Our results 
reveal that a very high concentration of group II elements in 
stanene destroys the linear dispersion of a bands (Dirac cones) 
at the Fermi level (Fig. 6 (c) and (d)). Therefore, our findings 
show that the best doping concentration of group II 1s 12.5%, 
because large energy gap was obtained at that percentage. 
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Fig. 7 Calculated Bands Energy of Stanene doped Meg: (a) 6.25% (b) 
12.5% (c) 18.75% and (d) 25% 


We also have computed and analyzed the total density of states 
(TDOS) and partial density of states (PDOS) of stanene doped 
Mg as can be seen in Fig. 8. The result of DOS helps to 
elaborate the nature of the band gap and PDOS provides 
details information about the origin of bands for valence and 
conduction bands. The s-orbital of Sn atoms contribute to 
lowest valence band while middle of valence band which 
occurred approximately within (0 eV to -5.0 eV) range is 
contributed from s-orbitals of Mg atoms and s-orbital of Sn 
atoms. In the unoccupied conduction bands above Fermi level 
(0.0 eV to 5.0 eV), the s-orbital of Mg atoms and p-orbitals of 
both Sn and Mg are the main contributors. In general, we 
noticed that p-orbital of Sn atoms are responsible for material 
properties near Fermi level. 


(a) Stanene doped Mg (b) Slanene doped Mg 
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Fig. 8 Calculated DOS (Density of state) and (b) PDOS (partial 
density of states) for the Mg atoms doped Stanene 
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IV. CONCLUSION 


In summary we have investigated the electronic and structural 
properties of stanene single layer doped Be and Mg by means 
of DFT. Our findings show that the best doping concentration 
is 12.5% for both Be and Mg because, large energy gap was 
obtained at that percentage. The electronic structure of Be 
doped stanene shows that the Fermi level is shifted towards 
the valance band edge when compared to pure stanene and this 
confirms that the Be doped stanene 1s p-type of material. The 
Dirac point of stanene locates at I’ shifted by 0.38 and 0.5leV 
for 6.25 and 12.5 % respectively, an energy band gaps of 0.27 
and 0.50 eV were obtained above the Fermi level for 6.25 and 
12.5% respectively, which indicate that the Be doped stanene 
could be a degenerate semiconductor. The band structure of 
Stanene doped Mg was also computed for 6.25, 12.5, 18.75 
and 25.0% of Mg doping 1.e.1 Mg atom, 2 Mg atoms, 3 Mg 
atoms and 4 Mg atoms is doped in the 64 atomic supercell cells 
of stanene. The band gap remains slightly above the Fermi- 
level and extents to 0.34 eV for 6.25 % and 0.43eV for 12.5 
%, in the case of 18.75and 25.0 % the Dirac cone disappear 
completely, an energy gap of 0.28 eV and 0.60 eV were 
obtained above the Fermi level for 6.25 and 12.5% Our results 
reveal that a very high concentration of alkaline earth elements 
in stanene destroys the linear dispersion of a bands (Dirac 
cones) at the Fermi level. Thus finally, doping and 
manipulation of such impurities make stanene a reachable/ 
useful material and remove its disability for its employment in 
optoelectronics applications. 
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